Enhancement of Loop induced 125GeV Higgs pair production 
through Large- Extra-Dimensions model at the LHC 

Sun Hao 1 : Zhou Ya-Jin 2 t 
1 School of Physics and Technology, University of Jinan, Jinan, Shandong 250022, P.R.China 
2 School of Physics, Shandong University, Jinan, Shandong 250100, P.R. China 



Abstract 

Based on the analysis of 5 fb 1 of data at the LHC, the ATLAS and CMS collaborations 
p> [ have presented evidence for a Higgs boson with a mass in the 125 GeV range. We consider the 

^ ■ 125 GeV neutral Higgs pair production process in the context of large-extra-dimensions(LED) 

model including the Kaluza-Klcin(KK) excited gravitons at the LHC. We take into account 
the LED effects coming from gluon-gluon fusion and quark-antiquark collision channels as 
well as their corresponding next-to-lcading ordcr(NLO) QCD loop induced corrections. We 
analyse their impacts on both the total cross section and some key distributions. Indeed, 
pp — > HH has the clear advantage of a lower standard model(SM) background compare 
to process like pp — > jj, though its SM prediction is very small, it is shown that the LED 
model raises the cross section of Higgs pair production compare to its SM prediction and 
enhance the transverse momentum(p^ ) and invariant mass(M^fij) distributions especially 
at high scales of p^ and Mhh- By including the NLO QCD loop corrections, the scale 
dependence of total cross section can be reduced obviously. Choose suitable decay modes 
like HH — > bb'jj or HH — s- and some simple cuts, we can strongly reduce the SM 

background but keep most of the LED effects, leading Higgs pair production a promising 
channel to search LED effects. 
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I. Introduction 



The hierarchy problem of the standard model (SM) strongly suggests new physics at TeV scale, 
and the idea that there exists extra dimensions (ED) which first proposed by Arkani-Hamed, 
Dimopoulos, and DvalipQ might provide a solution to this problem. They proposed a scenario in 
which the SM field is constrained to the common 3+1 space-time dimensions ("brane"), while 
gravity is free to propagate throughout a larger multidimensional space D = 5 + 4 ("bulk"). 
The picture of a massless graviton propagating in D dimensions is equal to the picture that 
numerous massive Kaluza-Klein (KK) gravitons propagate in 4 dimensions. The fundamental 
Planck scale Ms is related to the Plank mass scale M Pl = G N 1/2 = 1.22 x 10 19 GeV according to 
the formula M|j = 8ttM^ +2 R s , where R and 5 are the size and number of the extra dimensions, 
respectively. If R is large enough to make Ms on the order of the electroweak symmetry breaking 
scale (~ 1 TeV), the hierarchy problem will be naturally solved, so this extra dimension model 
is called the large extra dimension model (LED) or the ADD model. Postulating Ms to be 1 
TeV, we get R ~ 10 13 cm for 5 = 1, which is obviously ruled out since it would modify Newton's 
law of gravity at solar-system distances; and we get R ~ 1 mm for 5 = 2, which is also ruled out 
by the torsion-balance experiments [2J. When 5 > 3, where R < 1 nm, it is possible to detect 
graviton signal at high energy colliders. 

Both the ATLAS and CMS collaborations have reported a SM Higgs-like excess at around 
Mjj = 125 GeV. If a SM-like Higgs particle is discovered in this particular mass range, an 
important additional test of the SM electroweak symmetry breaking sector is the measurement 
of the Higgs self- interactions. At hadron colliders, the pair production of Higgs bosons plays a 
distinctive role in understanding the Higgs mechanism 0. As the triple self-coupling of Higgs 
particles is involved in such production thus provide the experimental reconstruction of the 
Higgs potential. Precise measurement of this coupling could therefore give more insight on 
the mechanism of electroweak symmetry breaking. Compared to that of a single Higgs boson 
production, the signal-to-background ratio could significantly improved. The invariant mass 
scale of the single Higgs production is fixed by the Higgs mass, of order only ~ 125 GeV. Thus 
their detection through heavy quark decay modes suffer from large QCD backgrounds. Further 
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more, the Higgs pair production can give various final states, depending on the decay modes 
of the Higgs boson. In the mass range of (120,130)GeV, the Higgs boson decay modes with 
the largest branching fractions are H — >■ bb and H — > W + W~ . The most probable decay mode 
for a pair of Higgs bosons is HH — > bbbb. However, this mode is challenging to search for due 
to the fact that it is difficult to trigger on, and that it competes against the QCD multi-jet 
backgrounds that possess overwhelmingly large cross sections. In general, QCD backgrounds 
can be suppressed with the existence of leptons and missing energy. In the SM frame, the 
promising channels are pp — > HH — > bbr + T~ (bbp~ p + , 6&77)[I] and pp — > HH —> bbW + W~^[ 
for Mh ~ 125 GeV. More detail analysis can be found later. 

Another important distinctive feature of the Higgs pair production at the LHC is that the 
effects of physics beyond the SM can remarkably enhance the cross section with respect to 
that of the SM. Phenomenological studies of Higgs pair production have thus been performed 
in the context of the fourth generation model [B], the littlest Higgs model [7] and the Universal 
Extra Dimensions model [5]. For the large extra dimensional models, the tree level diagrams 
mediated by the Kaluza-Klein gravitons lead to a large total cross section. Such new theoretical 
approaches have drawn extensive attention in ref[S] for a comparison between supersymmetry 
and LED models. Exchange of virtual KK graviton or emission of a real KK mode could give 
rise to interesting phenomenological signals at TeV scale [1U|, lllj . Virtual effects of KK modes 
could lead to the enhancement of the cross section of pair productions in processes, for example, 
di-lepton, di-gauge boson (77, ZZ, W + W~), dijet, ti pair[T2J Q3J US [03 [El [H] etc. Some of 
these calculations have also been expanded into NLO QCD loop level. A detailed calculation of 
Higgs pair production in LED model has been performed recent lv[18|. however, the QCD loop 
induced calculation based on the LED context is still missing. 

Experimentally, the CMS Collaboration has performed a lot of search for LED on different 
final states at ^/s = 7 TeV\19\ I21j . By combining the diphoton, dimuon and dielectron 
channels, lower limits are set on the effective Planck scale in the range of 2.3~3.8 TeV at 
the 95% confidence level|22j. These limits are the most restrictive bounds on virtual graviton 
exchange up to date. Based on the analysis of 5 fb" 1 of data at the LHC, the ATLAS [23] 
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and CMS [21] collaborations have presented evidence for a Higgs boson with a mass in the 125 
GeV range. We thus concentrate on the 125 GeV Higgs pair production related to the latest 
measurement with the effects of the LED models and find the characteristic distribution of it 
up to QCD loop induced level. 

This paper is organized as follows: in section 2 we present a brief introduction to the related 
theory, in section 3 we present the main contribution to the process and section 3 is arranged 
to present the numerical results of our studies. Finally we summarize the results in the last 
section. 

II. Related theories 

The LED model consists of the gravity sector and the SM sector. The manifold which gravity 
propagates, is not the ordinary 4 dimensional spacetime manifold, but R4xM, where M is a 
5-torus with radius R and volume Vs = (2irR) s without loss of physical significance. In our 
calculation we use the de Donder gauge. The relevant Feynman rules involving spin-2 KK 
graviton and the relevant vertices in the LED model can be found in Ref \W\ [TT| [25] . 
We denote the process as: 

PP ij Gkk — > HH + X(ij = qq = uu, dd, cc, ss, bb, gg) (2-1) 

where Gkk denotes the Kaluza-Klein(KK) gravitons. The couplings between gravitons and 
SM particles are proportional to a constant named gravitational coupling k = \/T6ttGn, which 
can be expressed in terms of the fundamental scale Ms and the size of the compactified space 
radius R by k 2 R s = 8-k(Att) s / 2 T(5/2)M s ^ +2 \ In practical experiments, the contributions of the 
different KK modes have to be summed up, so the propagator is proportional to i/(sij — mV), 
where = (pi + pj) 2 and is the mass of the KK state n. When the effects of all the KK 

states are taken together, the amplitude is proportional to ^ % = D(s). If S > 2 this 

ft Si J~ m a+ lt 

summation is formally divergent as becomes large. We assume that the distribution has a 
ultraviolet cutoff at A ~ Ms, where the underlying theory becomes manifest. Then D(s) can 
be expressed as: 
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The imaginary part I(A/y/s) is from the summation over the many non-resonant KK states and 
its expression can be found in Ref . [11] . Finally the KK graviton propagator after summing over 
the KK states is: 

®KK ~ D ( S "> {Vl^Vufi + VltfVva ~ Jj^VlJ-vVap) (2.3) 

By adding all the Feynman rules in the LED model and the propagator given above, we can 
get the amplitudes (.M) for the SM and virtual KK graviton exchange contributions as well as 
their interferences. The cross section integral for hadron-hadron collisions can be written as 

CT = E/ J^dxidx 2 G i/Pa (xi,ii})G 3/Pb (x 2 ,iM 2 f )- J avgfaJ^^l^^pS^ 1 d * n (2 ' 4) 



where Gij/ PatPb (x, /U 2 ) represent the gluon or (anti)quark parton density functions, p a ,Pb indicate 
(anti)proton, \ir is the factorization scale which can be chosen equal the renormalization scale 
fi r when the loop calculation is included. av g fac is the times of spin-average factor, color-average 
factor and identical particle factor. |.M n ,| 2 presents the squared n-particle matrix element and 
divided by the flux factor [2s(27r) 3 ™ _4 ]. denotes the n-particle phase space differential. The 
quantity s = x\X2S is the effective center-of-mass(c.m.s.) energy, and the sum Ylij runs m 
case of quark-antiquark incident partons over all possible quark-antiquark combinations (i,j = 
u,d,s,c,b,u,d,s,c,b). In case of gg initial state the sum has only one term with i,j = g. 
Typically the latest new parton distributions for collider physics CT10[26 has been used in our 
calculation. 

The n-body phase space differential d& n and its integral $ n depend only on s and particle 
masses m; due to Lorentz invar iance: 

$ n (s,mx,m2,...,m n ) = J d$ n (s, mi, m 2 , m n ) 

/n n 
S 4 ((P*+Pj) - X» II - ml)9{pl) (2.5) 

k=l k=l 

with i and j denoting the incident particles and k running over all outgoing particles(k = l,...,n). 
In our calculation we use BASES [27] to do the phase space integration and Kaleu[28] to cross 
check with each other. 
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Figure 1: Part of Feynman diagrams for pp — > HH in the SM (a, b, c) and the tree level quark- 
antiquark collision(d) and gluon-gluon fusion(e) diagrams in LED model, where q represents 
u-,d-,c-,s-, and b-quark while Gkk represents spin-2 KK graviton. 

III. Related Process 

We separate the total contributions into two parts. One is the pure SM effects and the other is 
the LED effects, simply denoted as a and a LED , respectively. 

In SM, there is gluon-gluon fusion channel contribute to SM predictions through top-quark 
loops, see Fig. QJa, b, c) for more details. There is also b-quark contribution, but it is small. 
Other diagrams include the change of the loop arrow in FigQJa, b, c) and cross change of the 
legs in Fig. QJc) which are similar thus not shown, totally ten Feynman diagrams contribute. 
Contributions from the quark-antiquark collision can be safely omitted in the light fermion mass 
limits except b-b fusion through t-channels. However, it is only less of 0.5 percent of gluon- 
gluon fusion contribution [6], and not considered here. Thus we can define this contribution as 

n SM = n SM LO op 
u u gg 

Now let's see the second part: contributions to the LED effects a LED . Several distinct 
contributions contribute to this part, include both the tree level and QCD loop induced level. 
We separate and highlight them into four different parts: 
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(b) 



(c) 




(d) (e) 

Figure 2: The QCD one-loop Feynman diagrams for the partonic process qq — > Gkk HH 
(a)-(d) and the counter term diagrams(e) correspond to Fig{T]^d), where q represents u-,d-,c-,s-, 
b- and t-quark while Gkk represents spin-2 KK graviton. 

• The tree level contributions defined as Gqq,gg BORN which include both the quark-antiquark 
collision and gluon-gluon fusion contributions, see FigQTd,e), directly connected by the 
effects of KK excitation of gravitons. There are two gravitational vertices thus such con- 
tributions are in the effective order 0(k 2 ). 

• The 0(a s ) level one- loop virtual corrections to the leading order process as well as the 
renormalization of the leading order cross section which we define as -^^ born ^ LEDloop . 
The QCD loop level Feynman diagrams of the quark-antiquark collision initial state part 
can be found in FigJ21 totally 4 loop diagrams and 1 counterterm diagram, while the gluon- 
gluon fusion contributions can be found in Figj3]with 13 loop diagrams and 1 counterterm 
diagram as well. These contributions are in the effective order 0(k 2 (x s ). 

• The SM one-loop prediction interference to the leading order process which defined as 
a LED BO RN®SM L oop. These contributions are in the order 0((y/Ka)a s ). We keep these 
contributions for a fully consideration and it may also interesting to see how large are 
these contributions. 

• Finally the real gluon and quark emission contributions. The Feynman diagrams can 
be found in FigHl The cross sections are defined as which are also in the order 



0(n 2 a s ). We denote these parts by: 



(1) q(pi)q(p2) ->■ G K k -> H(p 3 )H(p 4 )g(p 5 ) 

(2) q(pi)g(p2) -> Gaa -> H(p 3 )H{p 4 )q{p 5 ) 

(3) 9(pi)q(P2) ^ G KK ^ H(p 3 )H(p 4 )(q(p 5 ) (3.1) 



and 



(1) g(Pi)g(P2) -> Gaa -> H(p 3 )H(p 4 )g(p 5 ) 

(2) q(pi)g(p2) -> Gaa -> H{p 3 )H{p 4 )q{p 5 ) 

(3) g(pi)q(p 2 ) -> Gaa ->• H(p 3 )H(p 4 )(q(p 5 ) 



(3.2) 



refer to emission process related to quark-quark collision and gluon-gluon fusion sub-processes, 
respectively. 




Figure 3: The QCD one-loop Feynman diagrams for the partonic process — )• Gaa - ► irff 
(a)-(m) and the counter term diagrams(n) correspond to FigQJe), where q represents u-,d-,c-,s-, 
b- and t-quark while Gaa represents spin-2 KK graviton. 
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Figure 4: The tree-level Feynman diagrams for the real gluon/light-(anti)quark emission sub- 
process qq — > Gkk — > HHg related to Eq J3.2f 1) [Fig. (a)-(d)] . q{q)g — > Gkk — > HHq(q) 
related to Eq|3Jp,3)[Fig.(e)-(g)], gg -»• G KK -»• HHg related to Eqj3^1)[Fig.(h)-(k)], 
9Q(q) ~^ Gkk — > HHq((q)) related to Eq l3. 3( 2.3) [Fig. (1)] . where q represents u-,d-,c-,s-, b- 
and t-quark while Gkk represents spin-2 KK graviton. 

Finally, we write the total cross section as sum of the above definitions to the formula below: 



tot 



a 



„SM | -.LED 
a + a 



-SM L oop 



gg 

_i_ —LEDrqrn _i_ n .LED B oRN<2>LED L oop i —LEDrral 

' qq qq qq 



qq 



(3.3) 
(3.4) 
(3.5) 



, -LEDrqrn i rrLED BO RN®LED L oop _i_ ^LEDrorn^SMrqop _i_ rrLED REAL /q p\ 

+ "gg +°4a + °q<) + CT <7<7 * r.^n^J 



99 



n SM LOO P , -LED , LED 
a gg + a qqNLO + a ggNLO 

^SMlqop , , LED 



99 



99 



99 



+ ® ggLO + °ggLO LO + + AcT 99 



LED 



(3.7) 
(3.8) 



The reason to do this is to introduce the symbolic notations we use later and make it 
convenient for our numerical calculations. Eq 13.31 means we separate the total contributions into 
two parts. Eq 13.41 related to the SM contributions which include only UV and IR safe terms. 
Eq l3.5l and Eq 13.61 related to the quark-antiquark collision and gluon-gluon fusion contributions. 
Both are summed by two parts. One is the born level contributions the other are the loop 
induced LED contributions where we use Aa^®® to notate. We use cr^^f LO and <7g^p LO to keep 
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notation of Eg J3.5I and I3.6| respectively. Finally, we define the SM prediction plus the leading 
order LED contributions as a LO and the other terms are the shifted terms coming from the loop 
induced LED effects, see Eq l3.8l 

There exist UV and soft /collinear IR singularities in the calculations of a ®® BORN ® E LOOp . 
To remove the UV divergences, we need only the wave function renormahzation constants for 
the quark and gluon fields. We introduce the renormahzation constants 5Z^ L R for massless 
quark (q=u,d,c,s,b) fields and 5Za for the gluon field defined as 

^ q> L,R = (l + 5Z KLiR )h< ! ,L,R, Af = {l + 8Z A )*A% (3.9) 



In the modified minimal subtraction (MS) renormahzation scheme the renormahzation con- 
stants for the massless quarks are expressed as 

6Z, h L = -^C r (A uv - Am), SZ i , tJ , = -^C F (A uv -A, R ), (3.10) 

sz * = £(^-5"r T ") Aw+ S(^-5"f KTy ) A ™ < 3J1 > 

In the above equations \i T is the renormahzation scale, Cp = |, Ca = 3, Tp = i, nH v = 6 
corresponds to the six flavor quarks (u, d, c, s, t, b), whereas n* = 5 is the number of 
the massless quarks (u, d, s, c, b). Moreover, Ayy = r(l + euv){^Y uv and &ir = 
^j^F(l + ejji) (4:Tr) eiR refer to the UV and IR divergences, respectively. 

By adding renormahzation part to the virtual corrections, any ultraviolet (UV) singularities 
are regulated. Divergences arising from soft gluon emission removed by combining virtual and 
real emission corrections @qq t gg REAL - Singularities associated with initial state collinear gluon 



emission are absorbed into the definition of the parton distribution functions. We employ the MS 
scheme for the parton distributions functions. Similar to the virtual part, we utilize dimensional 
regularization to control the singularities of the radiative corrections, which are organized using 
the two cutoff phase space slicing(TCPSS) method [29"]. We adopt TCPSS to isolate the IR 
singularities by introducing two cutoff parameters 8 S and 5 C . An arbitrary small 5 S separates 
the three-body final state phase space into two regions: the soft region (E§ < 5 s \/ r §/2) and the 
hard region (E^, > 5 s ^/S/2). The 5 C separates hard region into the hard collinear (HC) region 



and hard noncollinear (HC) region. The criterion for separating the HC region is described 
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as follows: the region for real gluon/light-(anti) quark emission with I15 (or §25) < S c s (where 



§ij = (Pi + Pj) ) is called the HC region. Otherwise it is called the HC region. Then the cross 
section for each of the real emission partonic processes can be written as 

a = a + a = a + o + a . (,3.12) 

After combining all these contributions above, the UV and IR singularities in a tot are exactly 
canceled. Dependence on the arbitrary small cutoff parameters S s and S c are exactly vanished. 
These cancelations can be verified numerically in our numerical calculations. We use Feynarts 
to create the amplitudes and FormCalc to perform the numerical calculation. To compute the 
one-loop tensor integrals we use our modified LoopTools [30\ [3T| [32] which supply an automatic 
numerical checks for the Inferared(IR) divergence and cross check with OneLoop[33] at enough 
random phase space points. 

IV. Numerical results and discussions 

In the SM frame, given the small total cross section, it is clear that even for y/s = 14 TeV and a 
target luminosity of 0(1OOO/6 _1 ) one need to focus on the Higgs decay channels with the largest 
branching ratios to visible final states to observe pp — > HH + X such as H — > 66(59.48%), H — > 
W~W + {20.78%) , H r^r + ( / u-/^ + )(6.12%)[Ml- A feasibility study on the decay channels 
HH — > 6677, HH — > bbfi~ finding that with 600 fb~ 1 one expects 6 signal and 11 background 
events, giving s significance of about 1.5<r[l]. For the channel HH — > bbW~W + — > bblvjj, Ref[5] 
claimed with 57 signal and 119 background events at 600/6" 1 by employing new techniques and 
assuming good r reconstruction efficiency(~ 80%), thus make such channel a promising one. As 
can be seen, though challenge, Higgs pair production in SM can be measured. What's more, we 
hope to concentrate on the Higgs pair production under LED frame and find out whether this 
production can be enhanced and given more chance to be tested. 

In the numerical calculations, we take the input parameters as Mz = 91.1876 GeV, My/ = 
80.399 GeV[35], "(Mf) -1 = 1/127.934, M H = 125 GeV[23l[2l]. The factorization and renor- 
malization scales are chosen to be the same [i r = fif = /io = Mh = 125 GeV. Numerical analysis 
is done at the colliding c.m.s. energy 8TeV and 14 TeV LHC for the early and future LHC. 
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The UV and IR safeties are verified numerically after combining all the contributions at the 
QCD loop level. To check the UV and IR divergence cancelation, we display a random phase 
space point as well as the cancelation for different divergent parameters, see in Tabl^U One 
thing that should be emphasized is a LED BO R N ®LED LO op + a LED BO RN®SM LO op should include 
the conterterm contributions as well as the soft and collinear singularity terms coming from the 
hard emission. We implement this into our monte carlo codes which supply an automatic check 
of the dependence on these divergence parameters. We can see the UV and IR divergence can 
be canceled at high precision level in all the phase space, thus leading the continuance of our 
following calculation. 

A Random Phase Space Point for 

VP ~> KgQjCgg) H (P3)H{p 4 ), ij = (qq, gg) 

Pi=(152.42677337684606 152.42677337684606) 
p 2 =(152.42677337684606 -152.42677337684606) 
p 3 =(152.42677337684606 86.461267952196764 11.548609707441768) 
ff 4 =(152.42677337684606 -86.461267952196764 -11.548609707441768) 



divergence 


LEDborn 
L 11,99 


LED B ORN®LtiO LO op i ^LEDborn^SMlqop [fkl 
u qq.gg ' gg 11 


£uv ej R em ^ 

-±- = 4- = -±- = io 10 

e uv e rFf e/H 


0.8461630560055 
0.8461630560055 


0.1259656696 911 
0.1259656696 636 



Table 1: The UV and IR divergence cancelation at one given random phase space 
point for the loop contribution. Notice that in the loop terms, we include the 
conterterm contribution as well as the soft and collinear singularity terms coming 
from the hard emission contributions. 



Since the total cross section is independent of the soft cutoff 5 S (= AEg/Ef,, Ef, = \^§/2) and 
the collinear cutoff 5 C , we display the curves of the loop induced QCD corrections to integrated 
cross section for pp — > ij — > Gkk — > HH + X process versus the cutoff 5 S at the ^/s = 14 
TeV LHC in the LED model, where we take fif = fi r = fio = Mh, M s = 3.5 TeV ,5 = 3 
and 5 C = 5 S /100. For the strong coupling constant a s (fj,), we use the two-loop evolution of it 
with the QCD parameter A n f =5 = 226 MeV and get a s (fj,o) = 0.113. Nf is the number of the 
active flavors. Some of the LED results are listed in Table 2. The terms gg, uu, dd refer to pure 
gluon-gluon fusion, uu + cc, dd + ss + bb collisions respectively. We show the 5 S dependence on 
all the sub-contribution terms. It is shown clearly that the NLO QCD correction (A.a^^ NLO ) 
does not depend on the arbitrarily chosen values of d s and 8 C within the calculation errors. In 
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o s 






Ac7 NLo{P h \ 








UU 


dd 


99 


5 x 10" 


-3 


7.8433049751 13E-03 


3.524461919262E-03 


1.946637625053E-02 


1 x 10' 


-3 


7.846245116613E-03 


3.522936306130E-03 


1.943421303121E-02 


5 x 10- 


-4 


7.846382881667E-03 


3.528078248755E-03 


1.945856138503E-02 


1 x 10" 


-4 


7.846562713918E-03 


3.534350158736E-03 


1.942754069450E-02 


5 x 10' 


-5 


7.844722993687E-03 


3.532865728080E-03 


1.936325027843E-02 


1 x 10" 


-5 


7.793411078227E-03 


3.534430130933E-03 


1.933905866994E-02 



Table 2: The dependence of the loop induced QCD correction to the integrated 
cross section for the pp — > ij — > Gkk — > HH + X(ij = qq,gg) at the yfs = 14TeV 
LHC in the LED model, where we set \i r = fj,f = fio = Mh, M s = 3.5TeV, 5 = 3 
and S c = <5 S /100. The terms gg, uu, dd refer to pure gluon-gluon fusion, uu + cc, 
dd + ss + bb collisions respectively. 



the further numerical calculations, we fix 5 S = 10 -4 and 5 C = <5 S /100. To satisfy the unitary 
constraint, we adopt the cut \f§ < M s for the whole phase space. 

In Figf5j we show the scale(/i) dependence of the loop induced SM contribution (a SM ) : 
leading order born (<J LO ) and QCD loop corrected total (a tot ) cross sections as well as their 
seperate sub-contributions (c^^) in the SM and LED model at the yfs = 14 TeV LHC, and 
we define the corresponding K-factor as = a tot /a LO . There we take the input as M s = 3.5 
TeV and 5 = 3. From this figure, we can see: 

• For the SM prediction a SM , see the dot-dashed line in FigJSl though coming from the SM 
loop diagrams(Fig{T^a,b,c)), its cross section changes from 44 fb to 20 fb if the scale /x goes 
from 0.2 to 2 ^q, which implies obvious fi dependence. 

• For the pure LO and QCD loop induced LED effects, contributions from the gluon fu- 
sion {&g^NLo) ^ s a bout two times of it coming from quark-antiquark collision (c^^q). 
Separately, the crfjf/lo^ = I'll 99) reduce the [i dependence on the corresponding c,^fo 
contributions, see the dotted and dashed lines in Fig JS] for more details. 

• Sum all the contributions together, we get the fi dependence on the total cross section a LO 
and a tot with the definition in Eq 13.31 see t ne solid lines in Figj5j We can see the cross 
section value in the LED model varies from 95.4(65.9) fb to 48.2(50.9) fb for the LO and 
NLO contribution when \x goes from 0.2 /io to 2 \iq at the \fs = 14 TeV LHC. So that the K- 
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Figure 5: The scale(/i) dependence of the loop induced SM contribution (<r ) [dash-dotted line], 
leading order born (cr LO ) and QCD loop corrected total (a tot ) [solid line] cross sections as well 
as their seperate sub-contributions (o"^^) [dotted and dashed lines] in the SM and LED model 
at the *fs = 14 TeV LHC with fi = M H , M s = 3.5 TeV and 5 = 3. 
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factor in the LED model varies from 1.44871 to 0.946948. We see from this panel that the 
NLO QCD corrections in the LED model totally reduce the factorization/renormalization 
scale uncertainty obviously. In further calculations we fix fi = /jlq = Mh- 



PP-> HHat14TeV LHC 



-■■ \o=3 
5=5'"',,. N, 



3 a fl ft ft ft a .-. r-. B B B 



Figure 6: The cross sections and the K-factor(cr to '/cr i0 ) for the process pp — > ij — > Gkk 
HH + X(ij = qq,gg) in the SM(cj 5M ) and LED model at both the born(o- LO ) level and loop 
induced level(a to *) as functions of M s with fi = M H = 125 GeV at the y/1 = 14 TeV future 
LHC. The solid line refer to the SM prediction while the dotted and dashed line correspond to 
5 = 3 and 5 = 5, respectively. 



In upper panel of Figj6j we depict the cross sections for the process pp — > ij — > Gkk 
HH + X(ij = qq,gg) in the SM(a 5M ) and LED model at both the tree(cr LO ) level and NLO 
QCD loop induced level(cr to< ) as functions of the fundamental scale M s from 2.5 TeV to 6 TeV 
with /i = M H = 125 GeV at the y/s = 14 TeV future LHC. The solid line presents the SM 
prediction (a SM ) while the dotted and dashed line correspond to the total LED effects at both 
the born(<T LO ) level and loop induced level(<r toi ) with the extra dimension number 5 being 3 
and 5, respectively. Notice a tot = a SM + a LED and a LO are defined in Eq 13.31 and Eq 13.81 
respectively. From the figures one finds that the largest deviation from the SM due to LED 
occurs at small values of M s and 5. With the fixing value of 5, the born level LED effects and 
the QCD loop induced LED effects deviate at different value of M s behavior as the dotted and 
dashed lines presented. In order to compare the relative size of such deviation, we depict the 
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K-factor(cr to */V LO ) in the lower panel in Figj6l One can see for 5 = 3, the K-factor is larger than 
1 when M s < 3 TeV with the highest value 1.13, while M s > 3 TeV the K-factor can turn to 
it's lowest value 0.95 which means the loop effects reduce the born effects about 5%. For 6 = 5, 
the K-factor is always less than 1 if M s > 2.6 TeV. When M s > 4.5TeV, both line close to the 
the small value of 0.98. The same thing has been depicted in FigJ7]for the 8 TeV early LHC. 
We take the same input parameters as above, we found that though the LED effects enhance 
the SM prediction generally, the total loop induced effects reduce the born level LED effects. 
Their deviation become large if M s become small and the K-factor turn to their lowest value 
0.65 with M s = 2.5 TeV. When M s become larger then 5 TeV, no matter the value of 5, the 
K-factor turns to 1 which means the born level and loop induced LED effects become close to 
each other and also close to their SM predictions. 
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PP -> HH at 8TeV LHC 
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Figure 7: The cross sections and the K-factor(cr to */o" i0 ) f° r the process pp 



ij 



Gkk 



HH + X(ij = qq,gg) in the SM(cr 5M ) and LED model at both the born(<r LO ) level and loop 
induced level(cr to ') as functions of M s with /x = M H = 125 GeV at the *fs = 8 TeV early LHC. 
The solid line refer to the SM prediction while the dotted and dashed line correspond to 5 = 3 
and 5 = 5, respectively. 



In FigJHJ we show the transverse momentum(pr)[left panel] and the rapidity(y) [right panel] 
distributions of the Higgs boson for the process pp — > ij — > Gkk HH + X(ij = qq,gg) in 
the SM(da SM /dp%, da SM /dy H ) and LED model(da tot / dp% , da tot /dy H ) at the 14 TeV LHC. 
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Figure 8: The transverse momentum(px)[left panel] and Rapidity(y) [right panel] distribution of 
Higgs bosons for the process pp — > ij — > Gkk —> HH + X(ij = qq,gg) at 14 TeV LHC with 
M s = 3.5TeV and fi = fio = Mh- The solid line presents the SM prediction while the dashed, 
dotted, dash-dotted and dash-dash-dotted and dash-dash-dot-dotted lines refer to 5 = 2, 3, 4, 5, 6, 
respectively. 

There the results are for M s = 3.5 TeV at the fixed value 3 for the number of extra dimensions 
and obtained by taking the input parameters mentioned above. The LED effects gently raise 
the SM prediction at values of high p^ regions. Rapidity distribution is defined as ^ with 
rj = iln(Pi • q)/(P2 ■ q), where Pi and P2 are incoming proton momenta and q is the related 
Higgs boson 4-momenta. As we can see, the rapidity distributions in the LED model show 
significantly narrow peaks around y = 0, which implies the large contributions at high pt 
region. 




M HH [ GeV l M HH [ GeV l 

Figure 9: The invariance mass(M##) distribution of Higgs boson pair for the process pp —> 
ij — > Gkk — > HH + X{ij = qq,gg) at 8TeV[left panel] and 14 TeV[right panel] early and future 
LHC with M s = 3.5Te^ and /i = Mjj. The solid line presents the SM prediction while the 
dashed, dotted lines refer to 5 = 3, 5, respectively. The upper one refer to a tot while the lower 
one refer to a . 
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The distributions of the invariance mass(M#//) distribution of Higgs boson pair for pp — > 
ij — > Gkk — > HH + X{ij = qq,gg) has been displayed in Figj9]at the 8 TeV[left panel] and 14 
TeV [right panel] early and future LHC. Same conclusion can be obtained that the LED effects 
enhance the distribution at the high Mhh region since in these regions the LED effect dominant 
the total (SM+LED) distribution, as a result, more KK modes contribute with the increase of 
Mhh- The lower the dimension is, the higher value of Mhh deviation from its SM prediction as 
presented in Figj9j It will be also very interesting to compare the results at 8 TeV with that at 
14 TeV LHC. We can see that though 14 TeV LHC has much higher c.m.s. energy than 8 TeV 
LHC, the theoretical amplitude at higher c.m.s. LHC is hampered by the unitary constraint 
v§ < M s , thus the LED effect is much suppressed at 14TeV LHC than 8 TeV LHC. Besides, 
at the 8 TeV LHC, the total contribution a tot deviate from its leading order contribution a LO 
obviously. The deviations become much larger compare to the 14 TeV results. While at 14 TeV 
LHC, the deviations only appear in the middle of Mhh region between a tot and a LO when the 
dimension is fixed to be the same value. 

As have been told, the distribution shows the crutial character that the LED effects are 
mainly collected at high p^ and Mhh region while the SM backgrounds are leaving in the 
lower regions. We thus apply two strict cuts pr{H) > 500 GeV and Mhh > 1 TeV to strongly 
suppress the SM backgrounds while leaving the LED effects almost unchanged. We concentrate 
on two channels: pp — > HH — > 6677 and pp — > HH — > bb^C /U + . The main backgrounds 
come from pp — > ZZ(ZH)( with Z — > 66( / u~/_t + ) can be a background to h — > 66(/i _ /U + )) 
and pp — > tt — > bb^ jjL + vi>. By taking into account all these backgrounds and adopting the 
kinematic cuts mensioned above, we display the signal background ratio in Tabl^S] and 
Tabl^H corresponding to the two different decay channels with M$ and 5 as running parameter 
and luminosity equal (100, 600) [/6 _1 ] at 14 TeV future LHC. From these two tables we can see 
that Higgs pair production in LED model can give significant signal to background ratio when 
the scale M s < 5 TeV and dimension 5 < 5. By using the 6677 channel, one can test LED 
effects at M s = 5 TeV and 5 = 4 with > 3 while for bbfi~ channel with > 5. Both 
channels shows good character in testing LED effects. In the tables, ~ (<)1 means the -7= is 
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s 

Vb 

M s [TeV] 


<5 


6677 channel @ 
= 2 5 = 3 


14 TeV with £= 
5 = A 


(100, 600) 
5 = 5 


5 = 


6 


3.5 


26 


64 


12.4 30 


11 26 


8.2 


20 


6.5 


16 


4.0 


11.2 


27.4 


6.4 15.6 


4.4 11 


3.1 


7.5 


2.6 


6.3 


4.5 


5.2 


13 


2.8 6.8 


1.9 4.5 


1.4 


3.3 


1.1 


2.6 


5.0 


2.5 


6.1 


1.8 4.2 


1.3 3.2 


~ 1 


1.5 


~ 1 


1.2 


5.5 


1.3 


3.2 


~ 1 1.5 


< 1 1 


< 1 


~ 1 


< 1 


~ 1 



Table 3: for the cahnnel pp — > HH — > 6677 after taking into account all back- 
grounds and adopt the kinematic cuts. 



bbfT r fj+ channel ® 14 TeV with £=(100, 600) 



M s [TeV] 


<5 


= 2 


<5 


= 3 


5 


= 4 


5 


= 5 


5 


= 6 


3.5 


60 


146 


29 


70 


25 


61 


19 


46 


15 


37 


4.0 


26 


64 


15 


36 


10 


25 


8 


18.4 


6 


14.4 


4.5 


12 


29 


7 


16 


4.2 


10.3 


3.2 


8 


3 


6 


5.0 


6 


14 


4 


7.3 


2 


5 


1.4 


3.4 


1.1 


3 


5.5 


3 


7 


2 


4 


1 


2.2 


~ 1 


1.4 


~ 1 


1.2 



Table 4: for the cahnnel pp — > HH — > bbfi fi + after taking into account all 
backgrounds and adopt the kinematic cuts. 



close to(less than) 1, however, this can be improved if we apply the cuts more strict or wait for 
the luminosity to be higher at the LHC. It will be interesting that pp — > Gkk HH can also 
be a promising process in testing LED effects. 

V. Summary 

In this work, we present the a full treatment of the QCD loop induced corrections to the 125 
GeV neutral Higgs pair production process pp — > ij — > Gkk HH +X{ij = qq, gg) at the early 
(y/s = 8 TeV) and future (y/s = 14 TeV) LHC in the context of large-extra-dimensions(LED) 
model including the Kaluza-Klein(KK) excited gravitons. We investigate the dependence of 
the total corrected cross sections on the renormalization/factorization scale /x, and study the 
corrected distributions of the transverse momentary ), invariant mass(M##) and rapidity(y^) 
distributions. Our final results show that the LED model raises the cross section of Higgs 
pair production compare to its SM prediction and enhance the transverse momentum (p^) and 
distributions at high scales of pif and Mhh- By including the QCD loop corrections, the scale 
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dependence of total cross section can be reduced obviously. Choose decay modes like HH — > 6677 
or HH — > bbfj,~ /i + and some simple cuts, we can strongly reduce the SM background but keep 
most of the LED effects leading Higgs pair production a promising channel to search for LED 
effects. 
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